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Figure 1: A drawing, illustrating jet production in a proton anti-proton collision with the hard scattering 
process, initial state, final state radiation and hadronisation (jet fragmentation) including the underlying 
event. 

1. Introduction 

A remarkable property of quantum chromodynamics (QCD) is the property of asymptotic free- 
dom. It is reflected by the prediction of the renormalisation group equation (RGE) dependence of 
the strong coupling constant tts on the momentum scale. Experimental tests of asymptotic freedom 
require precise determinations of over a large range of momentum scales. Its determination 
discussed here [|]][^ is based on an inclusive jet cross section measurement with the D0 detec- 
tor [^. The inclusive jet cross section cp-o^s^x/ dpTd\y\ has been measured using the Run II iterative 
midpoint cone (R = 0.7) algorithm ^ in the energy scheme. 110 data points are measured as a 
function of the momentum transverse to the beam axis, pj in six equidistant regions of \y\ (absolute 
rapidity) spanning together a range of < \y\ < 2.4. 

The initial state of an inelastic pp collision is given by non-perturbative Parton Distribution 
Functions (PDF's) / which are determined empirically and follow a momentum scale evolution 
according to QCD. The partonic final state is determined by a convolution of perturbative QCD 
(pQCD) scattering amplitudes and the PDF's. The cross section for the hard scattering process of 
partons can be calculated via a perturbative expansion in orders of as 

^pert.(a,) = £ (g) /p(a,) ® /p(a,) . (1.1) 

n 

The pQCD calculations are corrected for non-perturbative effects of hadronisation and underlying 
event (see Fig. [T]). The prediction of the theory can then be expressed as 

^theory (O^s) = f'pert. (^i) " ^non-pert. ; 

(1.2) 

where the factor Cnon-pen. takes into account the non-perturbative corrections. At the same time mea- 
surements are corrected for detector effects. Comparisons between measurements and predictions 
of theory are then accomplished at the hadronic final state [Q]. 

Previous determinations of by means of hadronic jet cross section measurements have had 
large uncertainties in comparison to other determinations. They have not been very competitive but 
now the tools for substantial improvements are available. More precise theory calculations together 
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with more and better measured data lead to smaller uncertainties and competitive results. The jet 
energy scale has been continuously improved over many years and reaches at present a fractional 
uncertainty of 1 — 2% over a wide range of jet transverse momenta from 50 GeV to 600 GeV. 

Tevatron jet cross section measurements cover a wide kinematic range. In the kinematic plane 
of momentum fraction x and scale unique regions of phase space are covered, in particular when 
forward jets up to absolute values of pseudorapidity of three are included. There is also overlap 
with measurements from ep collisions. Therefore the D0 data allow precision tests of pQCD and 
an independent extraction of ag. 



2. Extraction of oc^ 

The used perturbative calculations are a sum of a full calculation to order &{(X^) [ next-to- 
leading order (NLO)] combined with ^{(X^) resummation of (2-loop) terms from threshold cor- 
rections [^]. Adding the 2-loop threshold corrections yields a significant reduction of renormali- 
sation and factorisation scale dependence which in turn leads to significantly smaller uncertainties 
in the theoretical calculations. These calculations have been ignored for a long time and became 
only popular and accessible since they are included into fastNLO [|n|] which provides fast recal- 
culations for arbitrary PDF's and is based on NLOJET-i-i- [12] [ pj| ] and code from authors of ref. 



[^]. The calculations are performed in the MS scheme ^ for five active quark flavours using the 
next-to-next-to-leading logarithmic (3-loop) solution to the RGE. The PDF's are taken from the 
MSTW2008 next-to-next-to-leading order (NNLO) parameterisations [^] [|l^]. The jet transverse 
momentum is taken for the renormalisation and factorisation scale. 

The determination of is accomplished from inclusive jet cross section data points by min- 
imising a function between data and theory. The as{pT) values are evolved via the 3-loop 
solution to the RGE to as{Mz). The central as{Mz) result is obtained by minimising with 
respect to a^ (Mz) and nuisance parameters for correlated uncertainties. 

The procedure described above to determine a, requires the knowledge of CJpert.(«s(^z)) as 
a continuous function of as{Mz) over a as{Mz) range which covers possible fit results and uncer- 
tainties. The MSTW2008 NNLO (NLO) parameterisations are available for 21 different values of 
a,(Mz) in a range of 0.107 - 0.127 (0.110 - 0.130) in steps of 0.001. CTEQ 6.6 parameterisations 
[g are available up to NLO for 5 different values of a,{Mz). The MSTW2008 PDF's with fine 
spacing in as{Mz) over a wide range are used together with cubic spline interpolations for the main 
result. 

The D0 inclusive jet production cross section measurement [^] entered already in the 
MSTW2008 PDF's. The correlations between experimental and PDF uncertainties are not doc- 
umented. Therefore they can not be taken into account in the extraction of and the cross section 
data points which provide strong constraints on the PDF's have to be avoided. Precision Teva- 
tron jet data dominate models for gluon densities at high momentum fraction x and start to affect 
gluon densities at x > 0.25. At leading order, di-jet events access x values of x„ = ^'"^ and 
Xfo = xt^—^-^—^ where xj is defined by xj = The mapping is less unique in inclusive jet 
production cross section data where the x- value is not fully constrained given a measured bin of pj 
and \y\. The full kinematics is unknown since it is integrated over the number of jets. The momen- 
tum fraction can be approximately constrained in treating the sub-leading jets as central {\y\ = 0). 
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The requirement x < 0. 15 removes all data points for which more than half the cross section is pro- 
duced at Xmax ^ 0.25. 22 out of the 110 data points fulfil this requirement (see Fig. ^ left). They 
are distributed in the kinematic phase space as follows: nine data points in the interval \y\ < 0.4, 
50 < < 145 GeV; seven data points in the interval 0.4 < \y\ < 0.8, 50 < pr < 120 GeV; four 
data points in the interval 0.8 < \y\ < 1.2, 50 < pr < 90 GeV and two data points in the interval , 
50 <pt< 70 GeV. 

The pQCD uncertainties due to uncalculated higher orders are estimated in varying the renor- 
malisation and factorisation scales a factor of two up and down. In the kinematic region used for 
the extraction of a, the two scales show positively correlated effects on the cross sections. As a 
conservative estimate both scales are varied into the same direction simultaneously. The resulting 
uncertainties can not be treated as Gaussian. Therefore the fit is repeated for the scale variations 
and the differences with respect to the central value of as{Mz) are added in quadrature to the other 
uncertainties. 



3. Results 

A combined fit to the 22 selected data points, grouped in nine transverse momentum intervals, 
yields as{Mz) = 0.\\6\^^^^^. The results are shown in Fig. § (right) as nine values in the range 
of 50 < pt < 145 GeV with their total uncertainties which are largely correlated among the data 
points. The results are complemented by DIS jet data from HERA. The as{pT) results obtained 
here are consistent with the energy dependence predicted by the RGE. The combined result for 



as{Mz) is consistent with the combination of HERA jet data [15] and the world average [16]. The 
largest contribution to the total uncertainty comes from the correlated experimental uncertainties 
which are dominated by jet energy calibration, jet transverse momentum resolution and integrated 
luminosity. 




P,{GeV) (GeV) 



Figure 2: Left: Inclusive double differential jet production cross section measured with the D0 detector 
The 22 selected data points for the extraction of ttj are indicated by the surrounding blue polygon. 
Right: The strong coupling constant as a function of pj (top) and at Mz (bottom). For comparison HERA 
DIS jet data are superposed. Error bars correspond to total uncertainties. 
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Vaiiation of the non-perturbative uncertainties by a factor of two up and down changes the 
central value by +0.0003 and -0.0010. It does not affect the uncertainty of the combined as{Mz) 
result. Replacing the MSTW2008 NNLO PDF's by the CTEQ 6.6 PDF's increases the central re- 
sult by only -1-0.5%, which is much less than the PDF uncertainty. Excluding the 2-loop threshold 
corrections and using pure NLO pQCD together with the MSTW2008 NLO PDF's and the 2-loop 
solution to the RGE yields as{Mz) = O.HOl^QQQjg. The small increase in the central value orig- 
inates from the missing ^{(xf) contributions which are compensated by a corresponding increase 
of Us- The difference to the central result is well within the scale uncertainty of the NLO result. 

4. Conclusions 

The strong coupling constant has been determined from the D0 inclusive jet cross section mea- 
surement ||3|] using theory prediction at NLO with the resummation of 2-loop terms from threshold 
corrections, i.e. at NNLO accuracy. The (Xs{pt) results are consistent with the predicted energy 
dependence of the renormalisation group equation. The combined result from 22 selected data 
points is as{Mz) = 0.116llQQQ4g. This is the most precise result at a hadron collider. 
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